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Temporary domain walls are produced during the bifurcation era of multi-stream inflation. The 
observational effects from such a domain wall to inflationary perturbations are calculated. A domain 
wall induces a preferred direction in the sky, affecting the angular distribution of perturbations. 
A correlated suppression of multipoles together with an alignment of multipole moments on the 
preferred direction are generated. Other observational aspects of multi-stream inflation, including 
hemispherical asymmetry and cold spot are also briefly reviewed. 



I. INTRODUCTION 



o 

(N 

t^ Recent CMB observations show nice agreements with the simplest inflation model, with nearly Gaussian and nearly 

O I scale invariant fluctuations. However, there are still a few anomalies, mainly as position space features, first discussed 

.^ based on the WMAP data (pL^ and references therein) and recently confirmed by Planck ^ at similar or improved 
confidence levels. Those anomalies include 

(N 

^^ • The angular power spectrum, when locally estimated at different positions on the CMB sphere, are different 

/^ EIH]. Especially the power spectrum on a hemisphere centered at co-latitude = 110° and longitude (j) = 237° 

f- \ is greater than that on the opposite hemisphere. ^ 

fH • Looking at more localized scales, a cold spot [6^, 7^ is found on the CMB map. The cold spot spans about 5° 

O , in the CMB map and is about 70/iK colder than average. Such a spot and the anomalous low temperature is 

' more than Act unlikely from pure statistical chance. 

^ •In terms of preferred direction for the fluctuation modes, an anomalous alignment of the quadrupole and 

^ octopole is also noticed [8, 9 . The directions for those multipoles coincide with each other up to an error 

■ ■ of 9°, corresponding to an 98% confidence level out of statistical fluctuation (or 13° and 96.7% CL from a 
different analysis). 

^ • The fluctuations at multiple moment / < 40 are smaller than expected by 5% to 10%, with a statistical signifi- 

ON cance of 2.5a to 3cr fTO^ ^ 

in 

o 



Those anomalies characterize some non-Gaussian features of the CMB temperature fluctuations. However, they do 
not seem to follow directly from nearly scale invariant non- Gaussian correlation functions in momentum space. Because 
■^sj" the estimators for non-Gaussianity with nearly scale invariance (see [11] for a review) are still very consistent with 
^^ Gaussian. For example, the observational bounds on the local, equilateral and orthogonal shaped non-Gaussianities 

are 12 , 



f]^l^^ = 2.7 ± 5.8, f^^f = -42 ± 75, /^^f ° = -25 ± 39 (68% CL statistical) . (1) 



Thus one expects the non-Gaussian features discussed above, if they exist and are not by pure chance, should be 
J^ either highly scale dependent, or manifest themselves better in position space. It is still unclear if some of those 
anomalies originate from an incomplete removal of galactic and/or extragalactic foregrounds. Also there is a debate if 
the anomalies are over-interpreted with too strong posteriori choices. Nevertheless those anomalies are likely to carry 
more primordial information than that in the SH initial [1 . Here we assume those anomalies comes from primordial 
physics and search for possible explanations from inflation. 

The inflation model we focus on is multi-stream inflation p!3HT7] . During multi-stream inflation, the inflation 
trajectory hit a bump on the way and bifurcates. Different inflationary local universes, and hence different patches of 
the CMB sky, follow different trajectories of multi-stream inflation. This situation is illustrated in Fig. [l] Considering 



^ One can alternatively ask if the CMB sky has a dipolar asymmetry in perturbation power spectra ^5, . 

^ The low / suppression manifests itself in momentum space. Also it is not necessarily related to non-Gaussianity as discussed below. 
Nevertheless let us list this anomaly together and we will come back to this anomaly in our calculations. 



that different trajectories are localized at different local universes before they return to the horizon, multi-stream 
inflation is efficient in producing position space features. 

In the remainder of the paper, we shall first review the known features from multi-stream inflation in Sec. [TTj and 
discuss new features from the temporary domain wall in Sec. Ill in the context of multi-stream inflation. 




FIG. 1: Multi-stream inflation. In this example, the inflation trajectory hit a sharp feature on a two-field potential and has 
to bifurcate. Different patches of universe follow different trajectories of inflation. The boundaries of those patches are domain 
walls during the bifurcation era. However, when the trajectories recombine, the tension of the domain walls, and thus the 
domain walls themselves vanish. 



II. MULTI-STREAM INFLATION 

The bifurcation behavior in multi-stream inflation can be modeled by different shapes of potentials. For example, 
the bifurcation can be realized by lowering the potential energy of the preferred trajectories 
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Or alternatively, one can put a Gaussian shaped spike to raise the potential energy of the forbidden region in between 
the preferred trajectories. This can be achieved by adding 



AV{4>, x) = M^ exp (-a^' - 6x') + mix ■ 



(3) 



onto a double field slow roll potential or quasi-single field potential to make the bifurcation. In both examples, a 
small perturbation is added to allow asymmetry of the two trajectories (the terms mfx)- 

Here we focus on model independent predictions of multi-stream inflation, and shall not discuss the features from 
departure of slow roll, or field self-interactions during turning of trajectory (see [18l |19] for those discussions based 
on potential (l2|)). The features of multi-stream inflation include ^: 



Small scale power asymmetry [TT: When the inflation trajectory bifurcates, the small scale power spectra at 
different patches of the sky are calculated by perturbation theory around different classical trajectories. Locally, 



^ Some of the features (especially the ones related to the cold spot) correspond to different parameter space of multi-stream inflation. 
Thus they may not be realized all together from one single bifurcation. 



the power spectrum can be still calculated as 
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(4) 



However, the local Hubble parameter and slow roll parameter e can be different along those trajectories. As a 
result, different patches of the sky admits different small scale power spectra. 

For example, when the bifurcation takes place on or before the first e-fold (corresponding to the largest scales) 
of observable inflation, the asymmetry in the power spectrum is hemispherical. This case is illustrated in Fig. [2J 

• Average temperature asymmetry [13]: The total e-folding number along different trajectories could also vary 
in multi-stream inflation. As a result, different trajectories admit different reheating time and thus different 
temperature. The difference in temperature follows from a combination of Sachs- Wolfe effect [20] and the 
^ A/"- formalism 



AT = -AAT 
5 



(5) 



To have this temperature difference under control, it is assumed that the total e-folding number along different 
trajectories are protected by an approximate symmetry, such that any e-folding number difference between them 
are highly suppressed. (Otherwise, the additional fluctuation can easily become too large.) This temperature 
difference is illustrated in Fig. [2] 
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FIG. 2: Asymmetries on the CMB induced by multi-stream inflation. The e-folding number difference between the blue and 
the orange patches brings an averaged temperature difference, whereas the difference in local expansion rate and slow roll 
parameter brings a difference in small scale power spectra between different patches. 

• Cold spot from last scattering surface [15] : A cold spot on the CMB may be produced by an extended version 
of multi- stream inflation^. Note that the cold spot is much smaller than the rest of the CMB sky. Thus the 
probability for the inflaton to take trajectory A (denoted by pa) or trajetory B (denoted by pb) has to be 
different. For example, pb ^Pa- 

Let the total e-folds of observational inflation be A/tot, and the bifurcation takes place at the e-fold A^bif (counted 
from the end of inflation) . Then the radius of the multi-stream bubbles is about exp(A/'bif — A^tot) x3000Mpc. The 
number density of the multi-stream bubble is pb exp(Atot — A'bif)- Again the density contrast is characterized 
by the e-folding number difference AA^ between different trajectories. To have a cold spot from the Sachs- Wolfe 
effect, an over density bubble and thus a smaller number of e-folds is needed inside the bubble. The detailed 
calculation of the bubble profile can be found in [15 . 

• Cold spot from inner void [15]: In this case, the dominate effects are the late time ISW effect [20 , the Rees-Sciama 
effect [27] and the Ostriker-Vishniac effect [28]. The qualitative analysis is similar to the above Sachs- Wolfe 
effect case, but now the multi-stream bubble is located in the large scale structure (for example z ^ 1) and a 
greater number of e-folds is required inside the multi- stream bubble compared with outside. 



For other mechanisms to produce the cold spot, see for example [24tl26] 




FIG. 3: Cold spot produced by multi-stream inflation. Either an overdensity region on the last scattering surface (a), or a 
void in the large scale structure (b), can be produced from a multi-stream bubble. 

• Constraining multi-field inflation on random potentials: In case the bifurcations are not under control, multi- 
stream inflation produce too large anisotropics. Thus those random bifurcations are tightly constrained by 
current observations. In the case of a two-field random potential, the bifurcation probability is estimated in 
[17], which is verified by order of magnitude using numerical simulations. The key observation is that on a 
random potential, the isocurvature direction has about half chance to be tachyonic. Along the inflationary 
trajectory, the mass term in the isocurvature direction can be written as a function of the inflaton field value: 

ml = ml{4>). (6) 

On a tachyonic potential the isocurvature fluctuations grow exponentially until the mass squared flips sign. 
There is a phase transition controlled by the averaged frequency of sign flip of m^(^). On the slower side 
of the sign flip, the isocurvature fluctuation has long enough time to grow and bifurcation becomes frequent. 
The bifurcation further enlarges isocurvature fluctuation, and convert the isocurvature fluctuation to curvature 
fluctuation by the e-folding number difference A^ = AA/" between different trajectories. 

On a two field potential, the constraint on the shape of the random potential is not very tight. However, one 
expects the constraint becomes much tighter (at least with linear growth of bifurcation probability) for a random 
potential of many fields. The scaling property of bifurcation probability as a function of number of fields is not 
yet calculated. 

• Temporary domain wall. During the bifurcation era (in Fig. [l] where the two trajectories are separate), the 
position space patches following different trajectories are separated by domain walls. If the domain wall is much 
thinner than the inflationary Hubble radius, the walls themselves do not significantly alter the inflationary 
background dynamics [16 . 

For an isocurvature field with potential V{x) = ~i^ffX^ + \^X^i the domain wall has thickness / ~ l/rriQ^ and 
tension a ~ rn^^/\. Thus we require TTieff ^ H near the bifurcation ^. For the toy potential, m^^ = M^ — g'^cj)^. 

The domain wall formation may be fast or slow, depending on the choice of parameters. The x fi^ld needs 
to roll a distance x ^ ^^eff/V^ to form the domain wall. The speed of motion, when tUq^ ^ H^ can be 
estimated from x ^ i^^m^e^^^^^. Thus a quick (i.e. shorter than one Hubble time) formation of domain wall 
corresponds to meff < — 2ProductLog(— A^/^/2) or meff > — 2ProductLog( — 1, — A^/^/2). Those solutions are 
valid for A < (2/e)^ c^ 0.29. The latter solution satisfies TTieff > H (thin wall). Here we have assumed the 
inflaton field (j) varies fast such that ?7ieff changes quickly enough into a symmetry breaking potential. The 
parameter regimes for fast and slow bifurcations are plotted in Fig. |5J 



^ If this is not satisfied, it remains interesting to see how the domain wall inflates together with the background for a period of time, and 
what observational features can be produced. 
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FIG. 4: Bifurcation probability on a two dimensional random potential |17j . 
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FIG. 5: Fast and slow bifurcations and domain wall formation process. Here we restrict our attention to A < {2/e)^ ^ 0.29 
(for A greater than this value, the analytical solution does not apply). The unshaded regime corresponds to fast bifurcation, 
where the domain wall is formed within one Hubble time. There are two fast bifurcation regimes. The upper one corresponds 
to a thin wall (compared with Hubble radius). The lower one corresponds to a thicker wall. 



For the time being, we restrict our attention to a fast bifurcation, with a thin wall, which corresponds to 

meff > -2ProductLog(-l, -A^/'^/2) . (7) 

In this case, the domain wall can be considered static during its life time. 

Also, those walls disappear when the different inflation trajectories eventually recombine, or reheat into the 
same thing, or has some rotational symmetry in field space ^. Thus those walls do not cause disasters for late 
time physics. 

However, it is still interesting to see if those temporary domain walls could bring new physics at the level of 



cosmic perturbations. One step towards this question is taken in Sec. Ill 



We thank Hassan Firouzjahi for pointing out the third possibility. 



III. PREFERRED MULTIPOLE DIRECTION 

During multi-stream inflation, a domain wall is temporarily generated between the two classical inflation trajectories. 
In this section, we consider inflationary perturbations with the presence of a planar domain wall. The validity of this 
section is not restricted to multi-stream inflation. Most of the calculation also applies for more general inflation 
scenarios with a thin domain wall. 

The domain wall solution can be understood in different frames. Here we consider two frames namely the domain- 
wall-comoving frame and the observer-comoving frame. In the domain- wall-comoving frame, with the domain wall 
located at z = 0, the stress tensor of the domain wall can be written as 



T/f = KMl diag(l, -1, -1, -1) + (j5{z) diag(l, -1, -1, 0) , 



(8) 



where A is the effective cosmological constant that drives infiation, and a is the tension of the domain wah. The 
Einstein equations with the above stress tensor can be solved by [23130] (see also [3T]) ^ 



ds-" = 



a'^{ri + p\z\)' 



{-drf + dz^ + dx^ + dy'^ 
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where 
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(r + i), /?: 
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Note that j3 takes value in the range — 1 < /3 < 0. This solution is illustrated in Fig. [6) that the acceleration in the 
z direction is reduced compared with a purely inflationary background. As we shall see, the reduced acceleration 
corresponds to a suppression of fluctuation along the z direction. 

This solution can be understood in Fig.[6J that in the z direction, which is perpendicular to the direction of the wall, 
the acceleration is reduced compared with that on the x-y plane. Note that this statement is observer dependent. In 
the observer-comoving frame, we get a different conclusion. Nevertheless, this domain-wall-comoving frame relates 
more closely to the observations. 




FIG. 6: Additional acceleration during inflation from a domain wall, in the domain-wall-comoving frame. 



^ Here we consider the planar solution. Domain wall solutions with alternative topology is known in [32| . which follows more closely from 



Alternatively, one can do a coordinate transformation to get to the observer- comoving frame. When z > 0, the 
metric can be written as [29l [30] . 

ds^ = ^ {-df + dz^ + dx^ + df) , (12) 

with the coordinate transformation 

fi = j{r]^f3z), z = j{z^f3r]), x = x, y = y, (13) 



where 7 = l/^/l^^P^. Thus P behaves as a Lorentz boost between the observer-comoving frame (12), and the 
domain- wall-comoving frame where the domain wall is static ([9]). In the coordinate (12), the position of the domain 
wall is i = P'yr]. Thus the domain wall is moving away from the observer. The metric in the z < part can be written 
down similarly, with a boost in the opposite direction /3 -^ —f3. 

Now let us consider how such a domain wall affects the fluctuations, which originate from sub-Hubble scales. We 



consider a spectator scalar field S(j) gravitationally coupled to the domain wall. Using coordinate (12), the calculation 
is standard, with 

d^k 



/d k 
-^e*^-^iifc(7))ak + h.c. , (14) 

(15) 




where 



and 

Mv) = ^^(1 + ikfi)e-'^^ . (16) 



Here the initial condition of the perturbation is chosen such that the perturbation locally minimize the energy density 
on the z > patch. This coincides with the Bunch-Davies vacuum on this patch. Also, note that the horizon crossing 
time is controlled by fj instead of r]. 

To relate to observations ^, we have to transform to frame ([9|. On super-Hubble scales {—kfj <C 1), the transfor- 
mation can be written as 

= I (03 e-- [ql + ql + ^Vs] - ^ ^a,(,) + h.c. , (17) 

where we have used that 7) :^ 0, and thus 77 c^ —j3z. Note that the equal 77 hypersurfaces should be used here because 
it is f] that controls the horizon crossing and frozen behavior . In the second line, we have changed the integration 
variable into q such that the different spatial directions have different weight: ki = qi, k2 = ^2 and ks = 7^3. Thus 
the Fourier space correlation function can be calculated as 

{S4>cS<P^') = {2iTfS\ci - q')^ (9? + ql + J^qly-' , (18) 

where we have used S^Ck — k.') = 7~^(^^(q — q'). Note that the result is symmetric under P -^ —j3. Thus this 
correlation function can also be extended into the z < regime. 

Clearly, the two point correlation function ([l8| has a preferred direction, perpendicular to the domain wall. Using 
the time delay formula (" = —H5(j)/(j)^ the power spectrum for (" is now ^ 



^ Here we assume the domain wall (before it disappears) is located crossing the center of our observable universe. Thus the domain wall 
rest frame is preferred by the symmetry of the system. 

^ Note that the time dependence of the inflaton background ^, together with the gravitational perturbations are not included in the 
calculation here. We expect the local calculation around horizon crossing not being sensitive to those complexities. On the other 
hand, a curvaton scenario [35H37 I can be considered to by pass the gravitational perturbations. In the case of a curvaton, the angular 
dependence of perturbations is not modified. Only an overall factor in front of P^ is changed because of different conversion mechanism 
from isocurvature to curvature. 



Note that there the Lorentz boost factor 7 > 1. Thus the probabihty to find the fluctuation along the qs direction is 
suppressed compared with that on the other two directions. This angular dependence is plotted in Fig. [7| 
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FIG. 7: The power spectrum with a domain wall, along different directions. The domain wall is located at the equator. Note 
that this plot shows the preferred direction for the fluctuation power spectrum, not the temperature fluctuation itself. The 
suppressed power spectrum (in red color) is more probable to show up. This can be related to the lack of power for the first a 
few multipoles. 



When averaged over the angular directions, the power spectrum is 
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Thus the angular-averaged power spectrum is not modified. Only that a preferred direction is generated. However, 
this angular average of the power spectrum is naive, due to the existence of a preferred direction. 

The power spectrum itself, even as a function of vector q, does not fully characterize the statistical properties of 
perturbations. Because the power spectrum only shows the size of fluctuation given one direction, but does not show 
a difference in probability to get fluctuations pointing to different directions. 

In the case with a domain wall, this preferred direction does exist. This is because the effect from a domain wall 
corresponds to a Lorentz boost, and a Lorentz boost transforms a sphere into an ellipse. In the observer-comoving 
frame, there is no preferred direction which has a greater probability to observe fluctuations. In the domain- wall- 
comoving frame, the /3 factor brings a preferred direction. To find the preferred direction and corresponding probability 
distribution, let us consider points on the x-z plane. Using polar coordinate in both frames, we have 



X = RsinO = Rsini 



z = R cos = ^R cos t 



(21) 



Again the mapping is on equal 77 hypersurface, because such hypersurfaces control horizon crossing scales. Here has 
probability density satisfying distribution 



Prob(l9) = -sin^ 



The probability distribution of is determined by 
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Prob(l9) = 
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(22) 



(23) 



The probability distributions are plotted in Fig. [S] Note that the 6 angle here is the direction along which the 
perturbation propagates. The orientations of the multiples (defined in m[9]) point to the perpendicular direction to 
^, which peak at directions on the domain wall. On the other hand, the (j) direction still has a uniform distribution, 
following from the symmetry of the configuration. 

It is interesting to note that the directions with close to and tt are more probable. Those preferred directions 
have suppressed power spectra, as plotted in Fig. [7[ Thus it is interesting to see if the preferred direction is related 



to a small-/ suppression of the power spectrum ^^. 

Based on the probability distribution function, one can further check the probability of mode alignments. The angle 
between two vectors (^i,(/)i) and {02-,4>2) is 



arccos [cos c/i cos c/2 + sm c/i sm 62 cos (02 



(24) 



The probability for two modes with direction vectors lying within 5 < 10°, as a function of 7 is plotted in Fig.[9J Note 
that here we do not map the direction vector to the two dimensional CMB sphere. Thus the probabilities are lower 
than the numbers listed in [2^. The figure shows an increase of probability with increase of 7. Note that if one takes 
large 7, the multipoles are suppressed correspondingly. It remains interesting to see if such a correlated suppression 
of low multipole moments and alignment better fit the CMB data. 
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FIG. 8: Probability distribution function of direction angle for different 7 factors. The 7=1 choice corresponds to the case 
without a domain wall. 
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FIG. 9: Probability of alignment for two multipoles. The probability is calculated numerically based on 10^ random samples 
per 7 value. Note that the probability of alignment becomes much higher for large 7. 



Before ending up this section, a few comments are in order here. 

• We have matched observational quantities with the quantities in the rest frame of the domain wall. The 
calculation is done via a boost from the frame (12). Thus there is a question if such a boost produce different 



Doppler dipoles on different sides of the domain wall. Such dipoles, if exist, induce huge perturbations and thus 



10 



In multi-stream inflation, in case that the inflation trajectories recombines at a later time, the higher I multipoles are no longer 
suppressed. 



10 

constrain /3 strongly. However, one should note that the physical clock field during inflation is not boosted by 
the domain wall (at least in the case of multi-stream inflation). Thus there is no anisotropy from the boost of 
the background. On the other hand, the perturbations crosses the horizon following a boosted clock. Nearly 
scale invariance of the perturbations protects the consistence of the above picture. 

• In [30 , the cosmological perturbations are also calculated. However, [30] focus on a different situation, in which 
the whole observable universe is on one side of the domain wall. Thus the observables considered in [30 are 
different from the ones considered here. As a related topic, the cosmological perturbation of the domain wall 
itself is considered in [34] . 

IV. CONCLUSION AND DISCUSSIONS 

To conclude, we consider how domain wall affects the amplitude and angular distribution of inflationary fluctuations. 
The domain wall generates a Lorentz boost. The boost does not modify the angular-averaged power spectrum of 
perturbations from a naive average. However, the fluctuations along the direction perpendicular to the domain wall 
is suppressed by 1/7^, and the fluctuations along the directions parallel to the domain wall is amplifled by 7. 

Due to the boost, the density fluctuation becomes more probable to propagate along the direction perpendicular 
to the wall (Fig. [8|. The preferred direction reduces the flne tuning of mode alignment (Fig. |9J which affects for 
example the alignment of the quadrupole direction with the octopole direction) . This alignment is correlated with a 
suppression in power spectrum. It is interesting to further investigate how much better those correlated effects flt the 
CMB data. 

In the original model of multi-stream inflation, the domain wall exists only temporarily. The different inflation 
trajectories eventually recombine and the domain wall decays during inflation. We expect that after the decay of the 
domain wall, the domain wall energy density gets diluted quickly. Thus we do not observe a wall-like over density 
regime on cosmological scales (the opposite possibility is discussed in [TT). 

It would be interesting to consider the true time dependence of the domain wall, which may correspond to a time 
dependent boost factor. Intuitively after the perturbations become classical and frozen on super-Hubble scales, the 
vanishing tension of the domain wall no longer modifles the frozen perturbations. On the other hand, the smaller scale 
modes which exits the horizon after the disappearance of the domain wall should remain standard and unmodifled. 
Only perturbations which exit the horizon during the bifurcation era are affected. We hope to address this issue in 
future work. 

In the calculation a spectator field is considered thus strictly speaking one still need to consider the conversion from 
this field fiuctuation to curvature perturbation. We expect the calculation can be directly applied to the infiaton, 
considering the perturbation of the infiaton should behave like a spectator field on sub-Hubble and Hubble crossing 
scales in the (^^-gauge (as in the case of the standard cosmic perturbation theory). This assumption needs to be 
verified. 

Also, only the planar class of domain wall solutions is considered here. It remains to verify if multi-stream infiation 
approximately generates the planar domain wall (which looks intuitive) , or other classes of domain wall solutions are 
more preferred. 

Finally, the essential mechanism discussed here is the boost generated by the domain wall. It is interesting to 
consider if there are mechanisms other than from a domain wall, which generate a similar boost and thus similar 
predictions. 
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